Topological Superconductors in Correlated Chern Insulators 
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In this paper, we realize a topological superconductor (TSC) in correlated topological insulator 
- the interacting spinful Haldane model. We consider the electrons on the Haldane model with 
on-site negative-U interaction and then study its properties by mean field theory and random- 
phase-approximation (RPA) approach. We found that in the intermediate interaction region, the 
ground state becomes a TSC with the Chern number ±2. We also study its edge states and the 
zero modes of the 7r-flux. 



I. INTRODUCTION 

Different from the traditional superconductor (SC), 
topological superconductor (TSC) always has the 
topologically-protected gapless Majorana edge states[l|. 
For TSCs, people cannot use the local order parameter 
to characterize them. Instead, it is the topological in- 
variable that plays the role to classify the topological 
properties of TSCs. According to the characterization 
of "ten- fold way" from random matrix [2l-j3l, there exist 
three types of TSCs in two dimensions: D-type chiral 
TSC, of which the topological invariable is the Chern 
number^, C-type chiral TSC, of which the topologi- 
cal invariable is also the Chern number and Dlll-type 
TSCQ, of which the topological invariable is Z2 topo- 
logical invariable. For the D-typc chiral TSC, a typical 
example is two dimensional Px ± ipy chiral p-wave su- 
perconductor. This TSC has exotic topological proper- 
ties, such as the topological protected chiral Majorana 
edge states and Majorana zero mode on 7r-flux[l[. In 
particular, a 7r-flux with Majorana mode becomes non- 
Abelian anyons due to its non-Abelian statistics [H. To 
realize a D-type TSC people always consider the system 
with s-wave pairing SC and strong spin-orbital coupling 
in strong Zeeman field For the C-type chiral TSC, 

a typical example is two dimensional dx ± idy chiral d- 
wave superconductor. This TSC also has the topological 
protected chiral edge states and zero mode on tt- 
However, till now people cannot realize the d^iidy chiral 
d-wave superconductor in a physical system. 

In this paper we will show another road to realize a 
TSC in an interacting topological insulator - the inter- 
acting spinful Haldane model. The Haldane model is a 
lattice model that illustrates integer quantum Hall effect 
without Landau levels Q. We consider the electrons on 
the Haldane's model with on-site negative-U interaction 
and then study its properties by mean field theory and 
RPA approach. We found that in the intermediate in- 
teraction region, the ground state becomes a TSC with 
the Chern number ±2. For the ground state there is a 
local order parameter and the elementary excitations are 
gapless phase fluctuations and gapped quasi-particle (an 



electron or a hole). To characterize its topological prop- 
erties, we study its edge states and the zero modes on the 
TT-flux. Thus we propose that the s-wave topological su- 
perconductor/superfluid may be realized by putting two- 
component (two pseudo-spins) interacting fermions on a 
honeycomb optical lattice. 

The paper is organized as below. In Sec. II, we start 
with the Hamiltonian of the interacting spinful Haldane 
model on the honeycomb lattice. In Sec. HI, We calcu- 
late the superconducting order parameter with mean field 
approach and get a global phase diagram at zero temper- 
ature. In Sec. IV, we discuss the topological properties of 
the TSC, including the zero modes of the tt-Aux and the 
edge states. In Sec. V, we calculate the phase stiffness of 
the TSC. In Sec. VI, we discuss the physical realization 
of the s-wave topological superconductor/superfluid in a 
honeycomb optical lattice. Finally, the conclusions are 
given in Sec. VII. 



II. THE SPINFULL HALDANE MODEL WITH 
ATTRACTIVE INTERACTION 

In this paper, we consider the Hamiltonian of the corre- 
lated Chern insulator - spinful Haldane model on hone 
comb lattice with the on-site attractive interaction as 

ini 



E 



(1) 



and 



where C/ > is the on-site interaction strength. 
fiil are the number operators of electrons with up-spin 
and that with down-spin, respectively, a are the spin- 
indices representing spin-up (ct =t) s-nd spin-down (a =], 
) for electrons, fi is the chemical potential and we only 
consider half-filling case by setting /i = 0. iJn is the 
Hamiltonian of the spinful Haldane model which is given 
by 
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Here t and t' are the nearest neighbor and the next near- 
est neighbor, respectively. To break time-reversal sym- 
metry, we introduce a complex phase into the second- 
neighbor hopping, and set the direction of the positive 
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phase is clockwise (\(j)ij\ = H' denotes an on-site 
staggered energy which is 



E 




FIG. 1: (Color online) The illustration of the honeycomb lat- 
tice 

Using the Fourier transformations, the electronic an- 
nihilation operators on the two sublattices are written 
into 



CieB,a 



= 2^e Ok., 

k 

1 \ " ik-Ri? 



(2) 



Ns denotes the number of unit cells. For free fermions 
(the on-site Coulomb repulsion U is zero), the spectrum 
is 



where 



and 



1=1 



(3) 



(4) 



ty 3 + 2 cos (VSky) + 4 cos (3fcx/2) cos {V3ky/2) 



Ik 



= 2t' ^ sin (k • bj 



there exist energy gaps near the points ki = ^^(1, -^) 



From the spectrum of free fermions, we can see that 
lere exist energy gaps h 

andka = -^(1, ^) as 

AE ^ 2e-6V3t' 
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When the energy gap closed, we can get the phase bound- 
ary. We can see that there exist two phases in the free 
electron case: the topological insulater with anomalous 
quantum Hall effect (we call it QAH) and the normal 
band insulator (NI) state. For each spin-component of 
fermions, the Chern number is defined by 



C = — [ (fk[n ■ d^n X dyii] = ±1 
47r Jfi 



(6) 



where n is defined as n = -py where 



d = {dx, dy, dz) 



and 



dx = Rc^k, 
dy = Im^k, 
4 = 7fe + £• 



Then at zero temperature (T = 0), due to spin rotation 
symmetry, we obtain C — ±2 in QAH and C = in NI. 
The cigcn-state of the two states are depicted in Figl2] 
and Fig|3l One can see that there is gapless edge states 
in QAH while no edge states in NI. 




FIG. 2: (Color online) The edge states of nornal insulator with 
open armchair boundary. The parameters are t' /t = 0.05, 
e/t = 0.3. 



The parameters ai , a2 and as are the nearest neighbors 
of the A sublattice (See Fig.l) which are defined as 

ai = ao(^,^),a2 = ao(i,-^),a3 =ao(-l,0) (5) 

and bi = a2 — a3, b2 = as — ai, = ai — a2. oq is the 
length of the hexagon side and is set to be unit. 



III. MEAN FIELD PHASE DIAGRAM AND 
TOPOLOGICAL QUANTUM PHASE 
TRANSITION 

When increasing the interaction strength, we get an 
s-wave pairing SC order. The SC order parameter of the 
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crgy in the reduced Brillouin zone as 



FIG. 3: (Color online) The edge states of topological insulator 
(or QAH) with open armchair boundary. The parameters are 
t'/t = 0.05, e/t = 0.01. 



s-wave SC is 



Then in the mean field approach, the Hamiltonian can 
be written as 



(7) 



In the momentum space, it becomes 



(8) 



where the basis vector ^'j^ is ( aj^.j. a_k^ b\^^ ^-ki ) a-nd 



-1k+e -UA -Ck 

u _ I -UA 7, - e ^k 

" -Ck -7fc - e -UA 

Ck -UA -7fc + £ 



(9) 



After diagonalization, we can obtain the spectrum of the 
quasi-particles 



Ek, = ±V(i7^l + V(t^^) + (10) 



and 



Ek, = ±y (|7fc| - \/{UAr + e^r + lap. (11) 

Then in mean field approach, we get the self- 
consistency equation to derive A by minimizing the en- 



2 I (EL 



i = — F 



•/(|7fe| + V(t^^)' + £')' + l^-P 



(12) 



(|7fc|-V(^A)^+£')' + iefcP 



To determine the phase diagram, there are four types 
of phase transitions : the quantum phase transition be- 
tween topological s-wave SC order with A 7^ and QAH 
with A = 0, the topological quantum phase transition 
between topological s-wave SC order with A 7^ and 
the normal s-wave SC state with A 7^ 0, the quantum 
phase transition between NI with A = and QAH with 
A = 0, the quantum phase transition between NI and 
the normal s-wave SC state with A 7^ 0. In particular, 
the topological quantum phase transition between TSC 
and normal s-wave SC is determined by the condition of 
zero quasi-particle's energy gap 



A/ = 2(3V3t' - y^(C/A)Ve2) = 0. 



(13) 



After determining the phase boundaries, we plot the 
phase diagram in FigU) In the phase diagram of FigHl 
there exist four different quantum phases : QAH, NI, 
C = ±2 topological s-wave SC and the normal s-wave 
SC. 
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FIG. 4; (Color online) The phase diagram of the case e/t = 
0.15 at T = 0. There exist four phases : the NI state, the QAH 
state, the topological s-wave state and the normal s-wave. 

From FigHl one can see that in the non-interacting 
limit (?7 = 0) , the ground state is a C = ±2 topolog- 
ical insulator with QAH for t' > 0.0288t and NI for 
t' < 0.0288t. At t' = 0.0288t, the electron energy gap 
closes at high symmetry points in the momentum space. 
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FIG. 5: (Color online) The dispersion of electrons for t' — 
0.0288f when U — 0. We can see clearly that in the high 
symmetry point the energy gap is zero and the dispersion has 
a Dirac cone. 
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FIG. 6: (Color online) The superconductor order A for the 
case of e/t — 0.15 and t' /t — 0.03. (■7)a is the critical point 
of the superconductor order phase transition. 



As a result, a third order topological quantum phase 
transition occurs between QAH and NT. See the disper- 
sion of electrons for t' = 0.0288t in Fig.5. With the 
increasing of the on-site Coulomb interaction strength, 
the ground state can be an s-wave superconductor. For 
t' > 0.0288i, the quantum phase transition between QAH 
and s-wave order is always first order which is denoted 
by the black line in FigHl 

In the region of 0.0288^ < t' < 0.035t, due to the 
jumping of the superconductor order, the QAH state will 
turn into normal s-wave state without gap closing. In 
Figini and Figl?! we plot the SC order parameter and the 
energy gap for the case of e/t = 0.15, t' /t = 0.03. While 
in the region of t' > 0.035i, with the increasing of the 
interaction strength, the QAH state will turn into the 
topological s-wave SC order after crossing a first order 
phase transition (black line in Fig|4]) and then turn into 
the normal s-wave order crossing a second order quantum 
phase transition (red line in Fig|3|). In Figl5]and FiglHl 
we also plot the SC order parameter and the energy gap 
for the case of e/t — 0.15, t'/t = 0.1. In the region 
of t' < 0.0288i, the quantum phase transition between 
NI and s-wave SC order is always second order which is 
denoted by the black line in Fig|4l With the increase of 
the interaction strength U, due to the smoothly change 
of the SC order parameter, the NI state will turn into the 
normal s-wave order after crossing a second order phase 
transition. In FigfTOl we plot the SC order parameter for 
the case of e/t = 0.15, t'/t = 0.01. 



IV. TOPOLOGICAL PROPERTIES OF TSC 
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FIG. 7: (Color online) The energy gap A/ for the case of 
e/t = 0.15 and t'/t — 0.03. (7-)a is the critical point of the 
superconductor order phase transition. 
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In this section we will study the topological proper- 
ties in the TSC state. The s-wave paired ground state 
exhibits a non-trivial topological property that can be 
signified by the Chern number C ~ ±2 in the momen- 
tum space. Because the pairing is s-wave, such TSC bc- 



FIG. 8: (Color online) The superconductor order A for the 
case of e/t — 0.15 and t'/t = 0.1. (y)a is the critical point 
of the superconductor order phase transition and (y)c is the 
topological quantum phase transition. 
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FIG. 9: (Color online) The energy gap A/ for the case of 
ejt = 0.15 and t'/t — 0.1. (7-)a is the critical point of the 
superconductor order phase transition and (y)c is the topo- 
logical quantum phase transition. 



FIG. 11: (Color online) The edge states of normal SC with 
"zigzag" open boundary for parameters t' jt — 0.1, e/t — 0.15, 
U/t = 3.5. 




FIG. 10: (Color online) The superconductor order A for the 
case of e/t = 0.15 and t' /t — 0.01. (y)a is the critical point 
of the superconductor order phase transition. 



FIG. 12: (Color online) The edge states of TSC with "zigzag" 
open boundary for parameters t' /t = 0.1, e/t = 0.15, U/t = 
3.125. 



longs to topological superconductor. In this section we 
calculate the edge states and the zero modes around the 
TT-flux. 



A. Edge states 

In this part, we study the edge states of TSC. The 
dispersion of the edge states of normal SC and that of 
TSC are plotted in Fig. 11 and Fig. 12, respectively. From 
them, we find there exist no gapless edge states in the 
normal SC. While for TSC there always exist two gapless 
edge states. 




FIG. 13: (Color online) The particle density around a 7r-flux 
in TSC. 
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B. Zero modes of tt-Hux 

After recognizing the properties of edge states, we turn 
to study TT-flux. A 7r-flux denotes half a flux quantum 
one plaquette of the square lattice, i^o = For 
a TSC on a honeycomb lattice, a tt-Aux on a plaquette 
is confined at zero temperature and cannot be real exci- 
tation. From the numerical calculations, we found that 
there exist two zero modes around each tt-Aux. And in 
Fig. 13, we plot the particle density around a tt-Aux in 
TSC phase. 

V. PHASE STIFFNESS 

In this section we study the phase fluctuations of 
TSC going beyond mean field method by random-phase- 
approximation (RPA) approach [l^. Now the SC order 
parameters become A{i) = |A|e*^'. For the fluctuated 
SC order parameters A(i), there are two bosonic modes : 
a gapped amplitude mode and the Goldstone mode (the 
mode describing phase fluctuations). 



In the imaginary-time path-integral representation 
(/3 = l/T, H ~ ks = 1), using the Hubbard-Stratonovich 
transformation, the partition function can be written 
as [12 

/[dc;dc<,][dA*dA]e-^('=-=-'^*''^) (14) 



rP I 

S= / dTy2{-A*A + c*^drc;+Ht + Ht' +H, 
Jo r 

- A(r,T)c<,(r,T)cs(r,T)- A*(r,T)c;(r,T)c;(r,T)}. 

After integrating over the electron fields, in the momen- 
tum space, the action can be written as 

^ = I E . « A;(g)A,(g) + 2/3^7V - ^ InG-^ 

where Arf{q) = A^^o + A^(g), A^.o = Aq and G^^ = 
Go ^ + with 



Go' 



f 7k 
-Ao 

-a 




-Ao 

+ 7k 




1UJ„ 



-a 



- 7k 
-Ao 



lUJr, 





-Ao 

- 7k 



+ sj 



and 







AA{q) 




-Asiq) 

-aU-q) 



(15) 



Here we have used 




ArM) = 


-i= f dxAJx)e-'">-'' 

Vn J 


A^i-q) = 




A„(x) 




a;{x) = 


^ 



Next, we investigate the Gaussian fluctuations around 
the saddle point. Using the expansion of the natural 
logarithm of 



TrlnG^i = ln(Go"^ + G^^) 



(16) 



TrlnGo 1 - J2 ^Tr(GoGr') 



we expand the action to second order of the fluctuation 
fields A^(g). This procedure leads to an effective action 
as 

S = So + S[ (17) 
where the effective action at the saddle point is 

^°^7?J2 + Wp^N - TrlnGo ' (18) 

and the effective action corresponding the quantum fluc- 
tuations reads 

S[ = ^TyY,iGoip)G^'{p,p- q)Goip- q)G^\p- q,p)) 



= 5l + 52 



where 



and 



1=1 



(19) 
(20) 

(21) 
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Here At = {A\{-q),AA{q),AB{q),A*g{-q)) denotes the 
fluctuation field and Q' {q) is given by the following ma- 
trix 



Q'{q) = 



where 



Q'lM Qriiq) Q'lM Q'lM 

Q'2i{q) Q'22iq) ^23(9) QM 

QM Q'M Q'M Q3M 

Q'iM QM QM QM 



Q'lM = ^J2^^MGoilip- q), 



Quiq) J^^Goi2{p)Goi2{p - q), 

p 

Qiaiq) = j^^Go?.2{p)GQii{p - q), 
p 

Q'liiq) = j^^GoA2{p)GQi:i{p - q), 
p 

Q'2i{q) = j^^Go2i{p)Go2i{p - q), 

p 

Q'22{q) = J^^Goii{p)Go22{p ~ q), 

p 

Q'2-i{q) j^^Gozi{p)Gq2a{p - q), 

p 

Q'2i{q) = ^5I'^o4i(p)Go23(p- g), 



(22) 



(23) 



(24) 



Q'iliq) j^^Gm:i{p)Goii(j) - q), 
p 

032(9) "^J^Yl Goi3{p)Go42{p - q), 
p 

Qsaiq) ^ l^Yl Gq33{p)Go44{p - q), 
p 

'934(9) = ^ X! Go43{p)Go43{P - 9), 



Qiiiq) = j^^Go24{p)Go3i{p " q), 

p 

'942(9) = J^Yl Gqi4{p)Go32{p - q), 
p 

'943(9) ^ l^Yl Go34{p)Gq34{p - 9), 

p 

(944(9) = ^ XI G'044(p)<^033(P - 9)- 



(25) 



(26) 



See the details of Go in appendix. 

For the term Qi j (9) j using the Matsubara summation 
formula as ^/I^J^ui ^/(i'^ni ^ E) = np{E), at zero tem- 
perature we may obtain 
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Qzjiq) ^ j^Yl Goki{p)Gamn{p - q) 

p 

1 1 Am , Bk 



Cm , Dm A„in 



iuJm + E- (k) ILOm ~ (k) iwm - iujn + i?+ (k - q) 



lUJm - i^n - £^+(k - q) \UJra " i^n + (k - q) 

iuj^ - iw„ - (k - q) 
1 AmB 

ran , V ^ 



AT -luj^ - £;+ (k - q) - S+ (k) ^ -icj„ - S_ (k - q) - S+ (k) 

1/ , 



-icj„ + £;+(k-q)H 


^E, 


-(k) 


CMBynn 






-iw„ - £;+(k - q) - 


- E_ 


-(k) 








-iw„ + £;+(k- q) H 


-E^ 


-(k) 



k 



-iw„ + 


E_(k-q)- 


Hi?+(k) 








-iw„ - 


i?_(k-q)- 


-i?_(k) 




DMCmn 




-iw„ + 


£;_(k-q)- 


t-£^-(k) 


iven by 



the Fermi distribution function which reads npiE) = 
l/(e^^ + 1). When the usual analytic continuation 
iLOn — > + iO+ is performed, and in static limit a; = 0, 
Q'ijiq) is reduced into 



1 i 



^ ^ 1 -i , , . 

^, , , _ -1 Afc;(k)i?„„,(k - q) + Efc;(k)A„„(k - q) ^/2 -i 1 I • ^^^^ 

y«.W-^[2^ iJ^(k-q)+iJ+(k) VO i 1 

_^ ^ Afc;(k )An«(k - q) + Bfe;(k)C,„„(k - q) 



E 



i?_(k-q) + i?+(k) 

CfcKk)S,„„(k-q) + i?fc,(k)A 
i?+(k-q) + i?_(k) 



Cfc;(k)Umn(k — q) + -Dfc;(k)Cm„(k — q) ^ a rotated basis, the matrix Q'{q) can be changed 

^ E-(k — q) + -E-(k) into Q{q). Now we have the action of the quantum fluc- 

(27) tuations as 



In order to obtain the phase stiffness, we can express 



A^(<Z) -X^(g)e''^"(«) - [X^{q)+i9r,{q)]/V2, where Si = $tQ(^)$ (29) 

(y5^(g), A,, (q) and are the real fields. A,, (g) and 6*^ (q) 2 ' 



essentially can be regarded as the amplitude field and the 
phase field, respectively. For the general cases, ipr^iq) is 



small[l3[. Then the field Ag in the new basis 



/ ^Aiq) \ 
OA{q) 

0B{q) 

\ >^Biq) I where = (A^(g), 6*^(9), 6'i3(q), As (<?) and Q{q) in this 
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rotated basis i 


■eads 






Qiiiq) 


= (Q'li + 


Q21 + 


Q12 + Q22)/2) 


QMq) 


= (iQii + 


1(321 - 


- iOi2 - i(322)/2, 


Ql3{q) 


= (-iQ'i3 


-iQ23 


+ iQ'l4 + iQ24)/2, 


Qiiiq) 


= {Q'l3 + 


Q23 + 


Q'm + Q24)/2, 


Q2l{q) 


= (-iQ'ii 


+ iQ21 


- iQ'l2 + iQ22)/2, 


Q22{q) 


= {Q'li - 


Q2I- 


Q'l2 + Q22)/2, 


Q2M 


= {-Q'l3 


+ Q23' 


f Q'l4 - Q24)/2, 


Q24{q) 


= (-iQ'i3 


+ 1^23 


- iQ'l4 + iQ24)/2, 


Q3i(<z) 


= (iQai - 


iQ4i ^ 


-iQ32-iQ42)/2, 


Q32(<?) 


= (-Q31 


+ Q41- 


^" ^32 ^ Q42)/2) 


Q33(<Z) 


~ ('333 ~ 


Q43 ^ 


Q34 + Q44)/2, 


Qaiiq) 


= (iQ33 - 


i<343 ^ 


- i(334 - i(344)/2, 


Qiiiq) 


= {Q3I + 


Q'ii + 




Q42{q) 


= (iQai + 


iQ4i- 


- iQ32 - iQ42)/2, 


Q-isiq) 


= (-iQ33 


-iQ43 


+ iQ34 + iQ44)/2, 


Qiiiq) 


= (Q33 + 




Q34 + Q44)/2. 



Next, we focus on the phase fluctuations. Upon inte- 
gration of the amphtude fields we obtain a phase-only 
effective action as 



5l - |5^0U(?)gp/.a.e(g)%), (31) 
9 



where 9^q) = {0A(q),dBiq))- The matrix Qphase corre- 
sponding to phase-phase fluctuation is given by 



QpHaseiq) - ( ) (32) 

V Wph2l Wph22 J 



where 



Qphll — ~{Q2lM[^Qi2 + Q24M21Q12 + (321^^12^42 + ^24-^^22*342) + Q22 
Qphl2 = "((32lM{iQi3 4- (324M21<3l3 + Q2iM[2Q43 + Q24M^2Q43) + Q23 
Qph21 = -((33lMnQl2 + Q3iM^lQi2 + Q3lM[2Qi2 + Q34M;^2Q42) + Q32 
Qph22 ~ — ((331-A^nQl3 + Q3iM2iQl3 + (331-^^12(343 + ^34^^22(343) + Q33 




Qiiig) -Qiijq) 

Z(q) Z(q} 

-Qii(q) Qiijq) 

Z{q) Z(q) 



(33) 



Z{q) = Qii(g)Q44('?) - Qi4(g)Q4i(<z). (34) 

Finally, the total phase-only effective action becomes 

s[^p/2Y,o\q)Q',nuse{qnq) (35) 

q 

where the matrix Qphaseil) reads as 



n' (n\ — I U ^ •^phll Qphl2 

^phase\q) ^1 r) 1 _L n 

V Wph2l JJ + Wph22 

(Q'phll Q'phl2 \ 
Qph21 Qph22 J 



(36) 



The zero temperature superfluid density (phase stiffness) 
PsiO) is obtained by examining the Goldstone mode ac- 
tion in the static limit. We numerically extract the zero 
temperature superfluid stiffness Ps(0) by identifying 



'phll^p/i22 



'3p/il2'3!: 



ph21 



[V3p,(0)/2]q2 (37) 



for |q| ^ [13. 

Now we derive the superfluid density (the phase stiff- 
ness) for the TSC at zero temperature Ps{0) = Q- Af- 
ter obtaining the phase stiffness Ps(0), the effective La- 
grangian of the phase fluctuations is given by 



Leff ~ ^Ps(0)(V^)' 



(38) 



where 9{x) denotes the phase fluctuations. See the results 
in Fig. 15. 
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3.15 



we get a laser-field-generated effective gauge field on this 
honeycomb optical lattice as that given in Ref . [T6| . 

In addition, we consider a strong interaction via Fes- 
hbach resonance technique_ff[/ = — J7 "'i-T"'«>^ where 

i 

U > is the on-site attractive interaction strength 
Thus, we may have an s-wave SF state by tuning 
the interaction between fermions via Feshbach resonance 
technique. 

Finally we get an interacting two-component fermions 
system in 2D honeycomb optical lattice of the Haldane 
model as 



(40) 



FIG. 14: (Color online) The phase stiffness of TSC order at 
T = 0. The green region is the TSC order. 



where fi denotes the chemical potential. 



VI. PHYSICAL REALIZATION 



VII. CONCLUSION 



In this end, we design an effective fermion model with 
C = ±2 topological superfluid (SF) as its ground state 
in an optical lattice. When two-component fermions 
with repulsive interaction are put into a honeycomb op- 
tical lattice, one can get an effective interacting Haldane 
model. It is easy to change the potential barrier by vary- 
ing the laser intensities to tune the Hamiltonian param- 
eters including the hopping strength (<-term), the stag- 
gered potential (e-term) and the particle interaction ([/- 
term). 

We first design an optical lattice of the Haldane model. 
The Haldane model had been proposed in the cold 
atoms with three blue detuned standing- wave lasers, 
of which the optical potential is given by V{x, y) = 
^j=i 2 3^ sin^[fci(a;cos0j + y sin^j) + 7r/2] where Vq is 
the potential amplitude, 9i = 7r/3, 62 = 2tt/3, 9^ = 0, 
and fci is the optical wave vector in XY plane [13 • When 
two-component fermions are put into this honeycomb op- 
tical lattice, we may get an effective Haldane model by 
applying the Raman laser beams 



e ^'^c 



We introduce a complex phase (t>ij {\4>ij \ = 7r/2) to the 
next nearest neighbor hopping, of which the positive 
phase is set to be clockwise. To design a complex phase 
of the next nearest neighbor hopping for two-component 
fermions generated by gauge field on the optical lattice, 
Raman laser beams in XY plane are applied with spacial- 
dependent Rabi frequencies as fio sv[v{kLX + ^)e'^'°'^ and 

r^o cos(A:_La; -I- f )e~''''^^ (fc_L = 27r/(3a)) where a denotes 
the length between nearest neighbor lattice sites. Then 



In this paper we studied the correlated Chern insula- 
tors by considering the electrons on the Haldane's model 
with on-site negative-U interaction. We obtained its 
properties by using the mean field theory and RPA ap- 
proach. We found that in the intermediate interaction 
region, the ground state becomes a TSC with the Chern 
number ±2. To characterize its topological properties, 
we studied its edge states and the zero modes on the 
TT-flux. In the end we gave a proposal to realize such 
s-wave topological superconductor/superfluid by putting 
two-component (two pseudo-spins) interacting fermions 
on a honeycomb optical lattice. 
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VIII. APPENDIX: THE MATRIX Go OF RPA 
APPROACH TO DERIVE THE PHASE 
STIFFNESS OF TSC 



In this appendix, we give the details about the matrix 

Go (= (Gq-I)-!) as 



Gn — 



/ Go, 11 

Go, 12 



G, 



0.13 



Go, 12 Go, 13 Go, 14 \ 
Go, 22 Go, 23 Go, 24 
Go. 23 Go, 33 Gi 



0.34 



(41) 



\ Go,i4 Gq 24 Go, 34 Go, 44 / 
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Each clement of Gq can be divided into where the term " numerator" is 

4 R 



numerator 



denominator 



(42) 



- ia;,„(Aij + By + dj + Aj)('-'m + + + + 7^) 
+ 2iw„,7k(^ij + By - - Dij)^Al +£2 

+ [(A,, - B,j)S+ + (Cj - Aj)^^-]Kn + lal' + A2 + e2 + ^2) 



-{A,, - B,,)E+ + (Cy - Aj)^-]27kV A§ + £2 



and the term "denominator" reads 

(iw„, + E+){\UJ,-n - E+){\UJm + E_)(lUJm. - E_). 

For the element Go,ij, we can express the coefficients 



as 



A, 



Cii — 





+ 




1 Qij 




4 




4B+ 


Ah, 


+ 




Pij — Qij 




4 




AE+ 


Mr, 






Pij + Qij 




4 




4B_ 


Ah, 










4 




AE_ 



(43) 



and the parameters for A/^, iVy , P^j, Qij are given as 
follows: for Gq.ii: we have 



= 1, 



Qii 

for Go 12, we have 



VATTi^' 

Pll = 7fc + 

-(£2 + A2 +£7fc) 



VA^ 



(44) 



12 



Afi2 = 



-An 



Pi: 



112 = 



v/A^T^ 
-Ao, 
Aq 7fc 

TaTT^ 



(45) 



for Go, 13, we have 



A/13 
Afi3 

^^13 
Ol3 



0, 

0, 



VAT 



for Go, 14, we have 



A^14 


= 0, 


iVi4 


= 0, 


Pli 


= 0, 


Qii 





a-Ao 



for Go, 22, we have 

A'/22 = 1, 
iV22 



for Go,; 



P22 = 7fe - £, 

'^22 — , .„ „ • 

-Go, 14, we have 

A/23 = 0, 
iV23 = 0, 
P23 = 0, 



Q23 = 



-gAo 
v/A^T 



(46) 



(47) 



(48) 



(49) 



12 



for Go, 24, we have 



M24 = 0, 
N2i = 0, 



for Go, 31 — Go.13 



M31 


= 0, 




= 0, 






Q31 





for Go, 32 = Gq 237 wc have 



Af32 = 0, 
N32 = 0, 
-P32 = 0, 



W32 = 

for Go, 33, we have 

M33 = 1 

^33 = 



-^Ao 



^33 = -(7fe + £), 



;33 



TAT 



for Go, 34, we have 



M: 



34 



0, 



Ao 



34 



VaTTI^ 

-P34 = -Ao , 

Ao 7fc 



^34 



VAf" 



for Go, 41 = Gq we have 



(50) 



M41 


= 0, 


Nil 


= 0, 


P41 


= 0, 







^Ao 



(51) 



for Go, 42 = Gg 24, we have 



(52) 



M42 = 0, 
7V42 = 0, 

A2 = 



Cue 



VAT 



for Go, 44, we have 



(53) 



A/44 = 1 

A^44 = 



— e 



VA^ 

^44 = -{ik - e), 

_2 



Q44 — 



A^ 



£7fc 



(54) 



(55) 



(56) 



(57) 



[1] G. E. Volovik, Zh. Eksp. Teor. Fiz. 94, 123 (1988) [Sov. 

Phys. JETP 67, 1804 (1988)]. 
[2] M. R. Zirnbauer, J. Math. Phys. 37, 4986 (1996). A. 

Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 

(1997). 

[3] A. Y. Kitaev, AIP Conf. Proc. 22, 1134 (2009). 
[4] Shinsei Ryu, Andreas P. Schnyder, Akira Furusaki, and 
Andreas W. W. Ludwig, New J. Phys. 12, 065010 (2010). 
[5] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 
[6] M. Sato, Yoshiro Takahashi and Satoshi Fujimoto, Phys. 



Rev. Lett. 103, 020401 (2009). 
[7] Chuan-wei Zhang, Sumanta Tewari, Roman M. Lutchyn, 

and S. Das Sarmal, Phys. Rev. Lett. 101, 160401 (2008). 
[8] F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 
[9] J. He, S. P. Kou, Y. Liang, S. P. Feng, Phys. Rev. B 83, 

205116 (2011). 

[10] J. He, Y. H. Zong, S. P. Kou, Y. Liang, S. P. Feng, Phys. 

Rev. B. 84, 035127 (2011). 
[11] J. He, Y. Liang, S. P. Kou, Phys. Rev. B 85, 205107 

(2012). 



13 



[12] E. Taylor, A. Griffin, N. Fukushima and Y. Ohashi, Phys. 

Rev. A 74, 063626 (2006). 
[13] M. Iskin and C. A. R. Sa de Melo, Phys. Rev. B 72, 

024512 (2005). 

[14] Erhai Zhao and Arun Paramekanti, Phy. Rev. Lett. 97, 
230404 (2006). 

[15] S. L. Zhu, B. Wang, L.-M. Duan, Phys. Rev. Lett. 98, 



260402 (2007). 

[16] L. B. Shao, Shi-Liang Zhu, Li. Sheng, D. Y. Xing, and 
Z. D. Wang, Phys. Rev. Lett. 101, 246810 (2008). 

[17] C. Chin, et al. Rev. Mod. Phys. 82, 1225-1286 (2010). 

[18] T. Kohler and K. Goral, Rev. Mod. Phys. 78, 1311-1361 
(2006). 



